Topological Phases of One-Dimensional Fermions: An Entanglement Point of View 



Department of Physics, 



Ari M. Turner 

University of California, Berkeley CA 94720, USA 



Frank Pollmann 

Department of Physics, University of California, Berkeley CA 94720, USA and 
Institute of Physics, Academia Sinica, Taipei 11529, Taiwan 



(N 

o 



OV 



u 
-t— > 



i 

o 
o 



> 
O 

o 



Erez Berg 

Department of Physics, Harvard University, Cambridge, MA 02138, USA 
(Dated: March 20, 2012) 

The effect of interactions on topological insulators and superconductors remains, to a large extent, 
an open problem. Here, we describe a framework for classifying phases of one-dimensional interacting 
fermions, focusing on spinless fermions with time-reversal symmetry and particle number parity 
conservation, using concepts of entanglement. In agreement with an example presented by Fidkowski 
et. al. (Ref. [if), we find that in the presence of interactions there are only eight distinct phases, 
which obey a Zs group structure. This is in contrast to the Z classification in the non-interacting 
case. Each of these eight phases is characterized by a unique set of bulk invariants, related to the 
transformation laws of its entanglement (Schmidt) eigenstates under symmetry operations, and has 
a characteristic degeneracy of its entanglement levels. If translational symmetry is present, the 
number of distinct phases increases to 16. 



I. INTRODUCTION 

Topological phases of matter are not characterized 
by a broken symmetry, but rather by an underlying 
topological structure that distinguishes them from other, 
topologically trivial phases. Such phases have attracted 
a great deal of attention recently, especially since the 
theoretical prediction^— and subsequent experimental 
observation^ of both two- and three-dimensional real- 
izations of time-reversal invariant topological insulators. 
These (as well as their predecessor, the integer quan- 
tum Hall effect) can be thought of as band insulators 
characterized by the topological structure of their Bloch 
bands. Similarly, topological superconductors^— are 
characterized by the topological nature of their fcrmionic 
quasi-particle spectrum. All these systems can be under- 
stood from a non-interacting point of view. A complete 
classification of all topological phases of non-interacting 
fermions, given their symmetries, has been given in Refs. 

MM 

In the presence of electron-electron interactions, the 
Hamiltonian cannot be reduced to a single particle ma- 
trix. Therefore, strictly speaking, the above classification 
scheme of topological phases cannot be used. Neverthe- 
less, in some classes of topological insulators, the topo- 
logical order has been argued to be robust even in the 
presence of interactions, by generalizing the correspond- 
ing topological invariant to the many-body caseii~— . In 
other classes, however, the situation in the interacting 
case remains unclear. 

In a recent breakthrough, Fidkowski and Kitaev stud- 
ied a one-dimensional model of spinless superconductors 
with time reversal symmetry— They found that in the 
presence of interactions, the free-fermion classification 
breaks down from Z to Zs, i.e. there are only eight 



distinct phases that survive in the presence of interac- 
tions (as opposed to an infinite number without interac- 
tions). To the best of our knowledge, this is the first case 
where the non-interacting picture in a class of topological 
phases is found to be radically modified by interactions. 
Ref. [l| constructs an explicit path in Hamiltonian space 
through which phases with different Z numbers mod(8) 
can be connected, and also discusses the stability of the 
edge states. However, a more general understanding of 
the classification of distinct phases in the presence of in- 
teractions (in particular, in terms of bulk properties of 
the ground state wavefunction) is left open. 

In this paper, we develop a framework for classifying 
phases of interacting fermions in one dimension based on 
bipartite entanglement of the ground states wave func- 
tion. The fact that entanglement is a useful quantity 
to probe topological properties of wave functions has 
been shown in several recent publications, see for exam- 
ple Refs. GUnl Our technique is based on a method 
which was introduced in Ref. [l8| for classifying phases in 
spin systems spins. This method has also been developed 
more fully and shown to give a complete ID classification 
by Ref. [l9| (at least when translational symmetry is not 
required). Here we generalize the method to fermionic 
systems. We find that the eight phases found in Ref. Q] 
are indeed topologically distinct, and characterize them 
in terms of a set of invariants. These phases cannot be 
continuously connected by any kind of interaction as long 
as time-reversal symmetry and fermion parity conserva- 
tion are preserved. 

The basic idea is to examine the behavior of the entan- 
glement (Schmidt) eigenstates of a segment in the bulk 
of the system under the symmetry group of the system. 
Topologically nontrivial phases can be recognized by the 
presence of "fractionalized" modes in the entanglement 
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spectrum, which transform differently under the symme- 
try group from the constituent microscopic degrees of 
freedom of the system (analogous to the half-integer spins 
at the ends of the spin one Heisenberg chain). The char- 
acter of the entanglement spectrum cannot change with- 
out a bulk phase transition, at which the nature of the 
ground state changes abruptly or the correlation length 
diverges. 

The behavior of the entanglement modes reflects the 
character of the physical topologically-protected modes 
at the boundary of the system. However, unlike the edge 
modes, the entanglement spectrum represent a truly bulk 
property of the ground state wavefunction, and as such, 
it is not sensitive to symmetry-breaking perturbations at 
the surface. 

We start in Sec. [TT] by introducing fermionic Hamil- 
tonians with pairing terms, through the example of a 
single Majorana chain model. The general framework to 
classify topological phases based on symmetry proper- 
ties of the entanglement eigenstates is presented in Sec- 
tion IIIII We apply it to fermionic systems with time- 
reversal invariance and fermion number parity conserva- 
tion, and derive the invariants characterizing the eight 
distinct phases and the degeneracies in their entangle- 
ment spectrum. These phases are shown in the next sec- 
tion to have a Zs group structure, defined through the 
rules for combining phases with different invariants. In 
Sec. IIV1 we demonstrate how to construct each phase by 
combining single chains. In Sec. [Vj we discuss the ad- 
ditional phases which arise if translational symmetry is 
imposed. The results are summarized and discussed in 
Sec. EH 

II. FERMIONIC MODELS WITH PAIRING 
TERMS 

We will investigate time-reversal invariant one- 
dimensional systems of spinless fermions, in which the 
particle number is conserved modulo two. (The classi- 
fication of topological phases is most interesting in this 
case.) Such a situation can be realized in a system in 
contact with a superconductor. As a simple example, 
consider the following Hamiltonian^: 

Ho = ^ ('VV : ' <7V : ' "•''•) ' " 'V V 

3 3 

(1) 

with t,u > 0. The operators c] {c-) create (annihilate) 
a spinless fermion on site j. The first term comprises 
hopping of fermions as well as the creation and annihila- 
tion of pairs of fermions while the second term acts as a 
chemical potential. The fermion parity operator 

Q = e™^i n > 

with n j = CjCj commutes with Hq as the total number of 
fermions -/Vtotai modulo two is conserved. Furthermore, 



the Hamiltonian is time-reversal symmetric. (For spin- 
less fermions, time reversal is represented by complex 
conjugation.) 

Let us begin by considering only the conservation of 
the fermion number parity. This symmetry of Hq allows 
us to distinguish two different phases. The system un- 
dergoes phase transitions at t = u, but no local order 
parameter can be used to distinguish the two phases on 
either side of the transition. However, they can be dis- 
tinguished by their edge states. In the phase u > t, the 
ground state for an open chain is unique while it is two- 
fold degenerate for t > «j2 If u = 0, one can check that 
these states are given by the equal weighted superposi- 
tion of all configurations with fixed fermion parity (i.e., 
an even number or an odd number of particles). The 
ground state degeneracy in this case can be understood 
in term of degrees of freedom at the two ends. The two 
ground states cannot be distinguished by any local ob- 
servable in the bulk, because in any finite region of either 
state, the parity can be either even or odd. However, the 
two states are distinguishable when the opposite ends are 
compared to one another: c' N C\ + cjycj + H.c. has a dif- 
ferent eigenvalue for the two states. Furthermore, we can 
transform the two ground states into each other by acting 
with either c\ +c\ or cn—c n on the two ends of the chain. 
In other words, there is a single fermionic state that is 
split between the ends of the chain, and the observable 
described above measures its occupation number. The 
two states are degenerate because the only distinction 
between them is long-range, while energy measures only 
local correlations. In the phase t < u, however, there is 
no such degeneracy. The picture provided remains true 
even if we include interactions as the arguments can be 
stated in a way that only requires the Hamiltonian to 
conserve the fermion parity. 

The edge properties have a simple explanation in a 
different representation defined by the following trans- 
formation: 

a, = Cj+c] (2) 
bj = -ifo-ct). (3) 

ctj and bj are Majorana operators; they obey the relations 
{a^aj} = {bi,bj} = <5y, {ai,bj} = 0, a i = aj and b t = 
bj. The fermion parity, 1 — 2rij of a site, is given by 
ibjaj. Using these operators, Hq can be written (up to a 
constant) as 

H ° = \ [} Yl h 3 a i+ 1 + U Y1 • ( 4 ) 

Observe that each unit cell now contains two operators. 
In the case where t = 0,u = 1, the ground state is de- 
scribed by iajbj = — 1, i.e., each site is vacant. In terms 
of these variables, the phase t = 1, u = is also simple: 
ibjaj+i = —1 when t is large. (This can be regarded 
as the parity of a fermion shared between sites j and 



3 




E S E 

FIG. 1: Illustration of a bipartition of a ID chain into a seg- 
ment (S) of length L and a surrounding environment (E). The 
operators Oa and Ob act on the edges of the segment. 

j + 1.) This requirement does not completely determine 
the ground state wave function in an open chain, though, 
because it leaves a\ and 6jv free. There are therefore two 
degenerate states characterized by the occupation of the 
fermion shared between the ends, ibpjai. 

The presence of time-reversal symmetry leads to ad- 
ditional distinctions between phases. Quadratic, time 
reversal invariant fcrmionic Hamiltonians with conserva- 
tion of the fermion number mod (2) have been shown^Jii 
to support phases classified by an integer n £ Z (class 
BDI, according to Ref. @). Each phase is character- 
ized by having n gapless Majorana modes at each edge, 
and the different phases cannot be smoothly connected 
to each other without closing the bulk gap. It was later 
founds that in the presence of interactions this classifi- 
cation breaks down to Zg. Now we will begin the main 
discussion, whose goal is to show how the eight distinct 
phases in the general (interacting) case can be under- 
stood and classified according to properties of their en- 
tanglement eigenstates under the symmetries of the sys- 
tem, namely time-reversal and fermion number parity 
conservation. 



III. CLASSIFYING PHASES BY SYMMETRY 
PROPERTIES OF THE ENTANGLEMENT 
EIGENSTATES 

In the preceding section, we discussed physical edge 
properties to characterize different phases. Below, we 
present an alternative method of classifying the phases, 
which involves only the bulk. This is achieved by exam- 
ining the "entanglement spectrum" 16 i 18 i 21 ~ ^H f a finite 
(but arbitrarily large) segmen t 14 ' 15 , embedded in an infi- 
nite system. The transformation law of the entanglement 
(or Schmidt) states under the symmetry group of the sys- 
tem can be used to distinguish between different phases, 
as we describe below. 

We consider a bipartition of a ID chain with periodic 
boundary conditions into a segment (S) of length L and a 
surrounding environment (E) of length N 3> L as shown 
in Figure [T] For the segment S, the reduced density 
matrix of the ground state wave function | ip) is given by 

P s = tr E mm- (5) 

It is convenient to define an "entanglement Hamiltonian" 



Hs such that 

Ps = e~ Hs . (6) 

Then, the low "energy" states of Hs are the most 
likely states of the segment S, when the entire system 
is in its ground state. We call the eigenstates |0 7 )s 
of Hs "entanglement eigenstates" and the eigenvalues 
E 1 "entanglement eigenvalues" . (These are the same as 

the Schmidt states defined by \ip) = e~~ l^^s^)^ 
where |X7)-E are the corresponding Schmidt states of the 
environment.) 

Our approach is based on an important observation for 
the entanglement Hamiltonian Hs- The low "entangle- 
ment energy" excitations of Hs in a d-dimensional sys- 
tem may be well described by a d — 1 dimensional effec- 
tive Hamiltonian (see also Ref. [HI for an interesting dis- 
cussion of this concept). We consider the entanglement 
spectrum of a sufficiently large one-dimensional segment 
S, and focus on the "low-lying" entanglement states with 
E~f < E cut (where -E cu t is an arbitrary constant). 

We now make a crucial observation: in a gapped sys- 
tem with a finite correlation length these states can 
only be distinguished from each other by their behav- 
ior within a certain distance of the ends of the segment 
S. This is justified by the following argument: suppose 
that we measure a correlation function C(£) = {OeOs)^, 
where Oe acts on sites in the environment E and Os 
acts on sites in the segment (far away from the edges), 
respectively. £ is the minimal distance between the sites 
on which Oe and Os act. We expect that as £ becomes 
large, C(£) — > {Oe)^{Os)iP- Let us take for Oe an op- 
erator that projects onto a particular Schmidt state of 
the environment, say \x~i)e- {Os can be any operator 
as long as it acts far from the edges of S.) When ap- 
plied to the ground state, Oe projects (through the en- 
tanglement between E and S) also onto the correspond- 
ing entanglement eigenstate \4>~ ( )s of S. Thus, we have 
(@s)<p y ~ (Cs) v- That is, (Os) in each eigenstate of Hs 
is the same as in the ground state; i.e., far enough from 
the edge of S, any eigenstate of Hs behaves essentially 
like the ground state. Therefore, an operator acting far 
from the edge of S cannot distinguish different Schmidt 
states; its expectation value must be the same in all of 
them. (A simple generalization of the above argument, 
using an off-diagonal Oe, shows that Os cannot connect 
different low-energy entanglement eigenstates if it is suf- 
ficiently far from the boundary.) 

Furthermore, one can show (see App. for the case 
of bosons) that any linear transformation applied to the 
subspacc of low-entanglement energy eigenstates of S can 
be represented by a local operator acting on sites within 
a distance £ from either boundary of S, with an accuracy 
that improves exponentially with £. In particular, every 
symmetry operation of the system is a unitary operation 
acting on the entanglement eigenstates,— and thus can 
be represented (within the low-entanglement basis) as a 
product of two operators Oa, Ob (Fig- U), acting on sites 
near the left and right edges of S, respectively. (Note 
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that Oa, Ob arc specific for a particular ground state 
Thus, Oa and Ob form a (projective) representa- 
tion of the symmetry group. Classifying the representa- 
tions formed by Oa, Ob distinguishes different phases, 
which cannot be adiabatically connected unless a phase 
boundary is crossed. At the phase boundary, either the 
character of the ground state changes discontinuously, or 
else the correlation length diverges and the two ends are 
no longer independent. 

Let us demonstrate these principles through the case 
of an SU(2)-symmetric integer spin chain. If the segment 
S is sufficiently long, there can be situations in which en- 
tanglement states come in degenerate setsi£ This can be 
seen from the fact that, according to the argument above, 
any SU(2) rotation can be represented accurately (within 
the low-entanglement energy subspace) in terms of two 
generators, and S^, which act only within a distance 
I of the left and right edges of S, respectively. Within this 
subspace, Stot ~ Sa + Sb, where Stot is the total spin. 
Both S^, Sb can be block-diagonalized into irreducible 
representations of SU(2) with well-defined angular mo- 
menta. Now, since Stot is an integer spin, there are two 
possibilities: either all the blocks in and Sb form inte- 
ger representations, or they all form half- integer represen- 
tations. This distinguishes two phases: a "trivial" phase 
in which Sa, Sb are both integer (e.g., a fully dimer- 
ized chain), and a "non-trivial" phase in which they are 
half-integer (such as the Haldane phase of the spin-one 
chair>22). In the latter phase, all the entanglement ener- 
gies of S must have a degeneracy of at least four—, due to 
the even degeneracy at each end, guaranteed by the pres- 
ence of the half- integer spin operators and S g . This 
is an alternative explanation of the symmetry-protection 
of the Haldane phase, discussed in Refs. [l8ll28l - l3ll . 



A. Fermion parity 

We now turn to discuss fermionic systems. Let us con- 
sider a Hamiltonian H that conserves the fermion par- 
ity Q as defined in ((TT| with Q 2 = 1. We show that 
one can distinguish two phases, a "trivial" phase and a 
"non-trivial" one. In the "non-trivial" phase, a segment's 
entanglement spectrum is doubly degenerate. The dou- 
ble degeneracy is related to a single fermionic degree of 
freedom, which is split between the opposite ends of the 
segment. 

Any eigenstate of H 7 and in particular the ground state 
is also an eigenstate of Q. Hence the resulting re- 
duced density matrix ps and the entanglement Hamil- 
tonian Jis both commute with Q. The eigenstates |</> 7 ) 
oil-is m ay therefore be classified by their Q-eigenvalues 
(q s = +1 if the fermion number Np is even and q s = — 1 
if N F is odd.) 

We now posit that it is possible to find an effective ex- 
pression for Q within the low-entanglcmcnt energy sub- 
space, of the form Q oc Q A Q B , where Q A and Q B are 
local operators which act near the left and right edges 



of S, respectively. This is analogous to the example of 
spin described in the previous section, in which the total 
spin can be represented as a sum of operators acting on 
the left and right edges. [In the fermionic case, the de- 
coupling of Q is multiplicative, since Q is itself a unitary 
symmetry rather than an SU(2) generator.] 

Now, Q A and Q B can have interesting relationships: 
Q A and Q B may be two fermionic operators (that is, 
each contains an odd number of creation or annihilation 
operators) or they may both be bosonic. Note that Q A 
or Q B cannot contain a sum of bosonic and fermionic 
terms, becasue Q oc Q A Q B is bosonic42 

If Q A , Q B contain an odd number of fermionic oper- 
ators, Q A Q B = —Q B Q A . To emphasize this "statistical 
correlation" between the two ends define an angle fj, = 
or 7r to distinguish between the two cases, so that 

Q A Q B = e^Q B Q A . (7) 

For (i = (tt), Q = Q A Q B (-iQ A Q B ). The factor of i 
must be introduced in the latter case for consistency with 
the anticommutation of Q A and Q B , given that Q 2 = 

(Q A ) 2 = (Q B ) 2 = 1. 

For example, consider the Hamiltonian in Eq. Q] with 
u = 0, t = 1. In this case, it is not difficult to show that 
there are only two entanglement eigenstates on the seg- 
ment S (with a nonzero weight in the density matrix), 
the same two states as the ground states of the origi- 
nal Hamiltonian restricted to this segment. The fermion 
parity of these states is given by 

Q = -iQ A Q B , (8) 

where 

Q A = ci + c\ 

Q B = -i(c L -c{). (9) 

We see immediately that Q A and Q B anti-commute, and 
therefore /it = 7r. Furthermore, [%s,Q A ] = [H-S,Q B ] = 
0. 

For any system with fi = 7r, the above commutation 
relations imply that all the eigenvalues of "Hs come in de- 
generate pairs. To see this, note that Q and Hs can be 
diagonalized simultaneously. Then, if %s\^>\) = E\\<j>\) 
and Q\4>\) = qx\<j>\) (where q\ = ±1), then the state 
IV'a) = Q^Va) is such that Hs\tp\) = E x \ipx) and 
Q\ip\) = —Q\\ip\), i-e. \ip\) is an independent eigenstate 
with eigenvalue E\. Indeed, for the Hamiltonian (Q]) with 
u = 0, we find a doubly degenerate entanglement level 
with E = log 2. (All other entanglement levels in that 
system have E = 00.) 

Note that, unlike the Haldane chain example in Sec. 
IIIIl the entanglement spectrum is two-fold (rather than 
four-fold) degenerate. This is a consequence of the fact 
that the degeneracy is not associated with either Q A or 
Q B alone; it is related to the occupation of the fermionic 
level formed by combining Q A + iQ B . i.e., it is shared 
between the two edges. 
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In a bosonic system, the states of the entanglement 
eigenstates can be represented by \4> 7 ) = \a(3). Here a, (3 
describe the states of the left and right ends of the chain; 
that is, they enumerate the eigenvalues of certain low- 
"energy" combinations of observables that are functions 
of I sites at the respective ends. This factorization is pos- 
sible for fcrmionic chains with /i = as well. However, for 
chains with n = tt, the extra qs variable describing the 
two-fold degeneracy cannot be written in terms of local 
observables belonging to either end (i.e., the fermionic 
degrees of freedom cannot be measured independently at 
the two ends of the segment). Therefore, the entangle- 
ment states should be labelled by three variables, \a/3qs), 
with the fermion parity qs of the entire chain represented 
explicitly.— 

In a nonintcracting system, the entanglement Hamilto- 
nian is also noninteracting . 32 ' 33 It can be represented in 
terms of entanglement modes. The only subtlety is that 
some of these modes may be Majorana modes, which 
satisfy m 2 = 1; rrv = in. There are two topologi- 
cally distinct phases, depending on whether there are 
an even or an odd number of Majorana modes at each 
end of the segment. The Q A and Q B operators de- 
fined above can be found explicitly. Given that the left 
edge has Nf low-energy fermionic entanglement modes, 
fA.a {pi = 1, ...,Nf), and N m Majorana modes, mA,p 
08 = 1, ... , N m ), Q A is given by 




and similarly for Q . Note that Q is a bosonic operator 
if N m is even (corresponding to fx = 0) and a Majorana 
operator if N m is odd (p = tt). 

When interactions are included, only the "total pari- 
ties" Q A , Q B are well-defined. The separate modes f a , 
mp can "decay" into other combinations of modes, but 
their total is closely related to the symmetry Q, as we 
have just outlined. 

B. Time Reversal Symmetry 

We shall now examine the consequences of time- 
reversal symmetry on the degeneracies of the entangle- 
ment energies. It turns out that the combination of time- 
reversal and fermion parity conservation can have non- 
trivial effects. 

Bosonic models. Let us introduce the approach by re- 
iterating the results for a bosonic chain in the presence 
of time reversal symmetry with [Hs,T] = 0. The eigen- 
states |</> 7 ) can be represented by |</> 7 ) = \a/3), where a, j3 
enumerate the low- "energy" states associated with the 
two edges. Entanglement eigenstates which differ in their 
a ((3) index can be connected by a local operator close to 
the left (right) edge, respectively. The transformation of 
the eigenstates of Us factors into parts referring to the 
two ends (see App. |A"|) . 



It can therefore be represented as a product of two 
unitary transformations U A ,U B acting on the ends of 
the segment, so that 

T\a,(3)=Y / U£ a U B f3 \a>p>) (11) 

a', /3' 

and [U A , U B ] = 0. For a discussion of subtleties related 
to the anti-unitarity of T, see App. [Bj Applying T to an 
eigenstate twice yields 

T 2 \a(3) = ]T U A a ,U B , p ,(U A y a , a (U B )l lf; \a"n- 

a',P',a",0" 

(12) 

Using T 2 = 1 and that the two ends of the segment are 
independent, it follows that 

U A {U A Y = U B (U B )* = exp(m)l, k = 0, tt (13) 

We can thus distinguish two different phases, correspond- 
ing to k = 0, tt. Let us now focus on the consequences 
for the entanglement spectrum. Assume that |</> 7 ) is 
an eigenstate of Hs with eigenvalue E 7 , then U A \4>j), 
U B \4>~f}, U A U B \(f>j) are also eigenstates with the same 
eigenvalue because Us commutes with Ua and Ub- If 
k = tt, the unitarics U A ,U B are anti-symmetric and 
thus the four states are mutually orthogonal, resulting 
in a four fold degeneracy of the entanglement spectrum. 
If k = 0, the entanglement spectrum does not necessar- 
ily have any degeneracies. For example, in the Haldane 
phase of spin-1 chains, we find k = it and therefore the 
entire entanglement spectrum of a segment S is four fold 
degenerate i 18 i 44 

This method may be generalized to give a classification 
of phases with any given set of symmetries. For each re- 
lationship between the physical symmetries (e.g., T 2 = 1 
in the case just described), there is a corresponding rela- 
tionship between the factored symmetries of the entan- 
glement spectrum 3 ^, in which certain phases (e.g., k) can 
appear. Certain combinations of these phases are "gauge 
invariant" (independent of how the phases of the fac- 
tored symmetries are chosen). These combinations dis- 
tinguish between topological phases. In fermionic mod- 
els, an additional possibility is that symmetry operators 
at opposite ends may either commute or anticommutc, 
as described in the previous section. 

Fermionic models. We now consider a Hamiltonian 
which has both fermion parity conservation with Q 2 = 1 
and time reversal symmetry with T = 1. In the pres- 
ence of both symmetries, we show that each of the two 
phases defined in the previous section (/i = 0, tt) can be 
subdivided into four different phases. Furthermore, we 
discuss the consequences for the entanglement spectrum 
in each case. As T simply takes the complex conjugate 
(spin degrees of freedom are not considered here), it does 
not change the total fermion number and thus [T, Q] = 0. 
We will now classify the phases by examining how the 
properties of the factored versions of Q and T depart 
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from the relations of the full transformations, T 2 = 1. 
[T,Q] = 0. 

We first consider the case [i = 0, i.e., [Q A ,Q B ] = 
0. Then, both Q A and Q B are bosonic operators, and 
can be diagonalizcd simultaneously. Then, we distinguish 
two cases: Q A T = e i<t 'TQ A with (f> = 0,tt, and similarly 
for Q B . Note that <j) has to be the same for Q A and 
Q B , since Q = Q A Q B satisfies [T,Q] = 0. If = tt, 
time reversal changes the parity of the fermion number in 
cither end. (A similar situation occurs at the vortex cores 
of time-reversal invariant topological superconductors^.) 
We now examine the two cases <fi = 0, tt separately. 

(/i = 0, <f> = 0, K = or 7r)- The case <j> = is anal- 
ogous to the bosonic case considered above, with two 
phases, one corresponding to k = tt, characterized by 
a four-fold degenerate entanglement spectrum of a seg- 
ment, and one to k — 0, in which there is no necessary 
degeneracy in the entanglement spectrum. 

(fj, = 0, <f) = tt, k = or tt)- If <p = it, U A 
and U B (defined through T = U A U B ; sec App. EJ) 
are both fcrmionic operators, since they change the 
fermion parity. We know that T 2 = U A U B (U A U B )* = 
-U A U A *U B U B * = 1. This can only be satisfied if 
U A (U A )* = exp(i«)l and U B (U B )* = -exp(i/c)l with 
k = 0,7T. Note that {Q A ,U A } = and {Q B ,U B ) = 0, 
where both U A and U B commute with Tis- {{■,■} de- 
notes an anti-commutator.) Therefore, each entangle- 
ment level is four-fold degenerate, where the degeneracy 
corresponds to states with all possible combinations of 
Q A = ±1 and Q B = ±1. 

Next, we consider the case \i = tt. In this case, 
{Q A ,Q B } = 0, so these two operators cannot be di- 
agonalized simultaneously. Rather, every entanglement 
eigenstate can be labelled by the eigenvalue q = ±1 of 
the parity operator Q = —iQ A Q B , where \a(3,q = ±1) 
are degenerate. Since [T, Q] = 0, we must have either 
[T,Q A ] = and {T,Q B } = 0, or vice versa. Therefore, 
we define a parameter (f> — 0, tt such that TQ A = e l ^Q A T 
and TQ B = e l ^ + ^Q B T. 

In this case (fi = tt), phases with <f> — and tt be- 
have very similarly. To see this, we just note that if 
4> = tt then QT commutes with Q A . Therefore, we will 
define a modification of time reversal that commutes with 
Q A , T' := QT if <j> = tt and T' := T if <j> = 0. Let 
the factors of T be T = U A 'U B '. One can check that 
U A , U B are bosonic. The entanglement spectrum can 
be divided into two sectors with a fixed value of q. The 
operator U A , being bosonic, depends only on a, (3 and 
hence acts the same way on both sectors. Define k by 

e in t = JjA'jjA'* = jjB'fjB'* The possible values for K 

are and tt: 

(/i = tt, <f> = or tt, K = 0)- If K = 0, each entangle- 
ment eigenstate in each q sector can be singly degenerate. 
Therefore, counting the q = ±1 degeneracy, each entan- 
glement eigenstate has a minimal degeneracy of 2. 

(/i = tt, (p = or tt, k — tt)- If k — tt, the spectrum 
in each of the ±g sectors is fourfold degenerate (for the 
same reasons as in the bosonic case above with k = tt, 



i.e., there is a Kramer's doublet at each edge). Taking 
the q = ±1 degeneracy into account, every entanglement 
eigenstate is at least eight fold degenerate. 

There are therefore eight different phases classified by 
triplets (fj,,4>,K), where each entry is or tt. [i = tt if 
there are Majorana modes at the ends of the segment, 
(j) = tt if time reversal and Q A anticommute, and k = tt 
if U A (U A )* = -1 (leading to Kramers' doublets at the 
edges of the segment). Since, as long as the time- reversal 
and fermion parity symmetries are preserved, (fx, </>, ft) 
can only take the values or tt, they cannot change 
smoothly; the only way for them to change is through 
a non-analytic change of the ground state wave function, 
i.e. a quantum phase transition. 

The eight different phases and their corresponding 
minimal degeneracies are summarized in Table U in Sec- 
tion [VI] The degeneracies illustrate a distinction between 
interacting and noninteracting systems. As we will show 
in Sec. IIV1 the eight distinct phases can be realized by 
taking M copies of a single chain in the large t phase, 
where M = 1,...,8. Without interactions, the degen- 
eracy of the Schmidt spectrum would be equal to 2 M . 
Interactions can partly lift this degeneracy but cannot 
connect the eight phases defined by (fi, <j), k) adiabati- 
cally. 

IV. ADDITION OF PHASES 

The eight phases we have just obtained obey a group 
structure, which is defined by the rules of "adding" them 
together. This group turns out to be Zs. The addition 
rules reveal interesting distinctions between bosons and 
fcrmions. Wc will work out the addition table in some 
detail. 

Two systems can be added together by placing them 
side by side: hence if one system is in phase Pi and an- 
other is in Pi, then the combined "ladder" system is in 
phase P\ + Pi- (The combined system then remains in 
P\ + P2 even when the two constituent systens are cou- 
pled, as long as the coupling Hamiltonian is symmetric 
under time-reversal and fermion parity, and the bulk gap 
does not collapse.) This rule creates a finite group for 
a given set of symmetries. In particular, every element 
in this group has an inverse and the trivial phase is the 
identity element. 

The distinction between fermionic and bosonic systems 
is related to the inverse operation. If the system consists 
only of bosons, then the inverse element of any phase P 
is its complex conjugate (i.e., its time reversal): 

P + P* = trivial phase. (14) 

For an example, consider a spin one Heisenberg chain. 
A single chain cannot be adiabatically connected to the 
trivial phase because its ends transform as spin 1/2 de- 
grees of freedom. However, as shown in Rcf. [3a|, two 
coupled chains can be connected continuously to the rung 
singlet phase (i.e., a product state of spin zeros on the 
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rungs). The two chains are no longer distinguished from 
the trivial phase by their ends because the two half- 
integer spins couple to form integer spin states. 

In general, phases of bosonic chains are distinguished 
by the projective representation of the symmetry groups 
acting on the entanglement eigenstates (see Ref. UM . 
Each element in the symmetry group ^2 is represented in 
the entanglement cigenbasis as a left-hand unitary matrix 
U A (E) acting on the left index of the state, and a right- 
hand matrix U B (E) acting on the right index. Then the 
combined operation of two elements Hi and E 2 is repre- 
sented by U A (ErEa) = e l P A ^^U A (Ei) U A (E 2 ), and 
similarly for U B . To see that P + P* is trivial, consider 
the eigenstates of the entanglement Hamiltonian for a 
segment in the combined system 

|a/3)cou P icd = p\a2^2) p* ■ (15) 

The left-hand matrix U representing a symmetry E for 
the coupled system is 

u^Xt ( s ) = ( s ) ( 16 ) 

where the second factor is complex-conjugated be- 
cause the second chain is time-reversed. Hence, the 
phase factors cancel, [7 A ' cou P led (Ei)J7 A > cou P led (E 2 ) = 
jjA,coupiod^vj 1 vj 2 ^ anc j the resulting system is in a trivial 
phase. 

Now we can try to build more complicated phases out 
of simpler ones by placing them sidc-by-side. For bosonic 
systems, this procedure does not generate new phases 
in the presence of time reversal symmetry. Time re- 
versal symmetry and Eq. p^|) imply that P + P = 0. 
Hence starting from one phase, it is not possible to get 
more than two phases (the original phase and the trivial 
phase) . There may be additional phases that would have 
to be built up from independent starting points. The 
group is always a product of ZVs, in other words. 

However, for fermionic spin chains, P and P* are not 
necessarily inverses. The p-wave superconducting state 
Pi described by Eq. ([1]) with t > /i which is an order 
eight phase, as discovered by Fidkowski and Kitaev-i, is 
an illustration. Eq. (fT6)) breaks down because operators 
on the two chains can anticommutc with each other. In 
fact, starting from a single Majorana chain, we can gen- 
erate all possible combinations of fi, </>, and n. We now 
demonstrate this idea for a number of examples: 

1. Consider the Majorana chain, the ground state of 
Eq. ([T} with t > u which is in the (/x, </>, n) = 
(71", 0, 0) phase. When two copies are combined to- 
gether, the resulting phase has (/i, (j), /i) CO upicd = 
(0,7r,0), i.e., the ends are not Majorana fermions 
any more, but time reversal changes the fermion 
parity of the ends. The fermion parity for the 
segment of the combined chain is given by Q = 
QiQ 2 = (iQ A Q B )(iQ A Q B )^ = (Q?Q£)(QfQ B ), 
where the Q A ,Q B are fermionic parity operators 
of a chain n with Majorana ends (see Sec. IIII Al) . 



One can measure the parities of the ends separately 
because Q A = —iQ A Q A and a similar operator on 
B are bosonic operators, so /i CO upicd = 0. On the 
other hand </> CO upied = tt because T anti-commutes 
with Q A on account of the factor of i. Furthermore, 
one finds that K C oupiod = t. (See App. [BT ) 

2. Consider two chains with (//, <p, k) = (0, 7r,7r). Q a , 
Q A are bosonic, therefore Q^ uplcd = Q1Q2 is 
bosonic as well, and //coupled = 0. Time reversal 
acting on the left edge is represented as U A nplcd = 
U A U A - It anti-commutes with Q A , Q A , but com- 
mutes with their product, therefore CO upicd = 0. 
Since both U A , U A change the fermion parity, they 
both have to be fermionic. Therefore {U A , U A } = 
and we get that 

^coupled ^coupled = iPl^^) iPl^^) 

= -ufu£*ufuf* 

= -1. (17) 

Hence ft C oupicd = 7T, and the resulting phase is la- 
belled by (0,0, tt). 

3. Combining two chains with (0, 0, 7r) finally gives the 
trivial phase, because all the symmetries are rep- 
resented by bosonic operators, therefore k simply 
doubles to give mod(27r). This conforms with 
the fact that the Majorana chain is an order eight 
element of the group. 

Working out the addition rule in general gives the table 
of phases which are summarized in Table [I] A concise 
way to describe the general addition rule is to define A = 
k + (f>(mod 27r). Then we represent a state by a 3-digit 

binary number ^, ^ J. These numbers add modulo 8 

when the phases are combined. 

V. TRANSLATIONAL IN VARI AN C E AND 9 

If, in addition to fermion parity conservation, transla- 
tional invariance is also present, the number of distinct 
phases is doubled. Below, we derive the associated invari- 
ant, 9, which can take the values or 7r, independent of 
the invariants (//, 0, n) described above. The degeneracy 
of the entanglement spectrum, however, is not modified 
in either the 9 — or ir phases, and is given by Table HI 

Let us consider a fermionic chain with translational in- 
variance. According to the Sec. IIII Al the fermion parity 
of a segment S = [1,L] extending from j = 1 to j = L 
(where L is much larger than the correlation length £) 
can be written as 

Q(l,L) = e-^f(L)Q A (l)Q B (L). (18) 

Here, we have kept track explicitly of the position of the 
operators Q A and Q B , and of an overall constant sign 
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/ (L) (which was absorbed into the definition of Q A and 
Q B before). Translational invariance removes the neces- 
sity of choosing the sign of Q B (and Q A ) separately for 
each segment. This symmetry also allows us to write the 
parity operator of the segment [1, L] in terms of those of 
the two segments Si = [1, L'\ and 5*2 = \L' + 1, L] (where 
L',L-L'->0 as 



Q(1,L) 



Q(1,L')Q{L' + 1,L) 
e-^f(L')f(L-L>) 



(1) 



(L')Q A (L' + l) 



(19) 
(L). 



Equating Eq. [TH] and QIJ] gives that, within the low- 
entanglement subspace, we must have 

/ (I/) / (L - L>) [e~^Q B (L') Q A (£' + 1)] = / (L) . 

(20) 

Eq. [20] can hold for every state in the low-entanglement 
subspace only if these states are all eigenstates of 
e-^Q B (£') Q A {V + 1). Q A and Q B can be defined in 
such a way that the corresponding eigenvalue is 1. Then, 
we get that / (L 1 ) f (L — L') = f (L), which is solved by 



f(L) 



A9L 



(21) 



From the requirement that [Q (1, L)] = 1, we get that 9 
can only take the values or it. Thus, each of the eight 
phases found in the previous section is further split into 
two distinct phases, corresponding to the two allowed 
values of 9. For example, for the Majorana chain model 
(Eq. [TJ, t = 0, u = +1 and t = 0, u = — 1 describe 
distinct phases, although both have [i = 0. The ground 
state has all sites occupied or unoccupied, corresponding 
to 9 = 7r or 9 = 0, respectively. 

Note that both 9 and /x have a concrete consequence 
not only for the entanglement spectrum but also for the 
parity of the ground state in periodic chains. If the length 
of the chain is much larger than the correlation length, 
the parity depends only on \x and 9 and the chain length, 
Thus, a phase with /i = 7r has an odd 
number of fcrmions on a chain of an even length.— The 
phase 9 determines whether the parity of the ground state 
alternates as a function of L or not. This is shown in 
appendix [Cj 



VI. SUMMARY AND OUTLOOK 

We have described a systematic procedure for classify- 
ing the phases of ID interacting fcrmions. We focussed on 
spinless fcrmions with time-reversal symmetry and par- 
ticle number parity conservation. In the non-interacting 
case, these models are classified by an integer number, 
i.e., by We used concepts of entanglement to clas- 

sify the phases in the presence of interactions. We de- 
rive an effective description of the dominant entangle- 
ment states which then allows us to recognize "topolog- 
ical" features based on projective representations of the 



symmetries. We found, in agreement with the results 
of Fidkowski et. al. (Ref. [lj), that in the presence of 
interactions there are only eight distinct phases. Each 
of these eight phases is characterized by a unique set of 
bulk invariants (fj,, <p, k), which can take the values or tt. 
These invariants are related to the transformation laws 
of the entanglement eigenstates under symmetry oper- 
ations, and the phases have a characteristic degeneracy 
of entanglement levels. Furthermore, the phases obey a 
Z§ group structure and each of the eight phases can be 
generated by adding single chains together. All possible 
phases and the addition rules are summarized in Table [I] 
If translational symmetry is also present, the number of 
distinct phases increases to 16. 

The symmetries we have focused on describe only one 
of the 10 Altland-Zirnbauer classes^! of topological insu- 
lators. The framework described here can also be used 
to show how the phases in these other classes are mod- 
ified by interactions. To analyze each of the classes of 
topological insulators (and in both interpretations when 
particle-hole symmetry is present), one only has to de- 
termine the appropriate algebra of symmetries and then 
determine the possible projective representations of this 
algebra. 

Interactions cause the meaning of the Altland- 
Zirnbauer classes to bifurcate, however. At the mean- 
field level, a superconductor has an emergent particle- 
hole symmetry in its band structure. Thus, the classes 
which have such a symmetry can be interpreted as de- 
scribing either superconductors or systems that have a 
true particle-hole symmetry (such as the Hubbard model 
for fermions with spin on a bipartite lattice at half- 
filling). When interactions are included, these two inter- 
pretations are distinct. Thus, the class BDI, for exam- 
ple, has particle-hole and time reversal symmetry. This 
can be interpreted as describing superconductors. This 
means that one symmetry, particle conservation, breaks 
down, and only fermion parity Q is left. The only two 
symmetries are T and Q, giving the problem treated here. 
BDI has an alternative interpretation, according to which 
it describes systems with a true particle-hole symmetry 
C that reverses the sign of (n*) — n. (Here, rii is the oc- 
cupation number of a site and n is the mean occupation 
number.) In this case, particle number N is conserved, 
and besides this there are two other symmetries T and C . 
These satisfy the algebra T 2 = C 2 = 1, CN + NC = 2nL 
where L is the length of the system. (Every other pair of 
these symmetries commute.) The set of phases is differ- 
ent for the two interpretations; particle number conser- 
vation rules out Majorana Fermions. A complete classifi- 
cation of systems in all Altland-Zirnbauer classes in one 
dimension, following cither of the two interpretations of 
particle-hole symmetry mentioned above, would be an 
interesting project for future work. 

A generalization of these results to higher dimensional 
systems is an interesting (and challenging) open prob- 
lem. In some of the symmetry classes of topological 
insulators and superconductors, strong arguments have 
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Number of chains 


(H, cj>, k) 


Degeneracy of segment 


1 


Or, 0,0) 


2 


2 


(0,7r,7r) 


4 


3 


(7r,7r,7r) 


8 


4 


(0,0, tt) 


4 


5 


(7r,0,7r) 


8 


6 


(0,7T,0) 


4 


7 


(7T,7T,0) 


2 


8 


(0,0,0) 


1 



TABLE I: Degeneracies and addition table. All possible 
phases of fermions are realized by simply taking copies of 
some number of Majorana chains (see next section); the first 
column is the number of chains. The next column gives the 
parameters classifying a given state. The third column gives 
the degeneracy of the Schmidt spectrum, counting both ends. 

been given that the non-interacting classification does 
not change when interactions are included. This is par- 
ticularly clear when the topological invariant is related to 
a quantized physical response, e.g., in the integer quan- 
tum Hall effec t 12 ' 13 and in 3D time-reversal invariant 
topological insulator*^. However, for other classes, the 
situation is less clear. For example, the non-interacting 
classification of 3D chiral superconductors is Z^J^ simi- 
lar to the one-dimensional case considered here. It would 
be interesting to consider the effect of interactions on the 
phase diagrams of such systems. 

Note As we were writing this article, we learned that 
a similar classification is being worked out by Fidkowski 
and Kitaevi^ Our results are consistent with theirs. 
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Appendix A: Low Energy Operators and the Ends of 
the Chain 

An intuitive argument, given above, suggests that low- 
energy operators acting on the entanglement cigenstates 
may be represented approximately by operators located 
near the ends of the chain: in each Schmidt state \ot(3)s 
the expectation values of the spin and other operators 
have some particular spatial dependence near the ends of 
the chain depending on a and (3, but this decays exponen- 
tially to the ground state away from the ends. Therefore 
it should be possible to transform between these states by 
using operators defined on just the ends. A special case 
is the effective representations of symmetries in terms of 
operators at the ends of the segments, which we used to 



define the topological phases. 

To see that these effective operators exist, one can use 
a matrix product state representation^ of the wave func- 
tion. (We focus here just on bosonic systems. For sys- 
tems including fermions a similar argument can be devel- 
oped using bosonization but the discussion of this general 
result gets complicated by the classification of the phases 
by f- ) 

A basis for the low-energy operators can be constructed 
as follows. For each fixed choice of a\ and ct2 let 
A {ct2 1 OL\) be the operator that transforms a\ into cti. 

We will now give an approximate representation for 
A (ai,a 2 ) that gives the correct transformation of low 
energy states oiUs- Let the matrices T mi A define the 
bulk state (m varies over a basis for the physical Hilbcrt 
space) . The ground state wavefunction of a ring of length 
N is given by 

\tp) = ^2 tr(r mi Ar„ l2 . . . Ar mjv A)|mim 2 • ■ • m N ) 

. (A1) 

T m , A can be brought into a canonical form, satisfying 
£ m r m Art n = £ m rt„Ar m = 1, where A is a diagonal 
matrix with non-negative entries ! 40 ' 41 For a generic wave 
function, T m , A are infinite-dimensional. 

Let us define the states \a(3)\L of a segment of the 
chain stretching from 1 to L, where L < N, as 

\afi)i,L = ^ ( r ™i Ar '"2 • • ■ Ar mj -) Q(3 |mim 2 . ..m L ) 
{">.} 

(A2) 

When L is large these states are nearly orthonormal, 

that is \(a'/3'\afi)i,L — 8 a 'a5p'p\ ~ Ce" , where £ is the 
length scale for the decay and C is a constant dependend- 
ing on the indices. In this limit, |a/3)ii are the entangle- 
ment eigenstates of the segment. On the other hand, if 
the length of the chain is fixed and a,f3,a',f3' increase, 
the orthonormality must eventually break down for high 
enough a, (3, a', (3' (since the Hilbert space of a segment 
of length L is finite). Indeed, C grows as a function of 
a, a', (3,(3'. 

The ground state wavefunction (|A1[) can now be writ- 
ten as 

af3) L +i,N- (A3) 

a,fi 

This gives the Schmidt decomposition into states of the 
environment \a(3) l+i,n and states of the chain \ap}i s L, 
with a Schmidt eigenvalue A Q A^ [or, equivalently, an en- 
tanglement energy E = E a + Ep = — 2(lnA Q + lnA^)]. 
The Schmidt eigenstates become orthogonal to each 
other in the limit N — > oo and L — > oo. 

We can now give an "effective" expression for O a in 
terms of local operators acting on sites !,...,£ (near 
the left edge of the segment valid for a low 

entanglement-energy subspace with E a < E cut . E cut is a 
cutoff which depends on I. The accuracy of our effective 
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expression improves as £ becomes larger (provided that 
N > £). Define 

= E MiAaiTlh* (A4) 

7=1 

where x is a cutoff of the entanglement spectrum that 
satisfies E 1= i x < E cut . O a ^ acts only on the £ first 
sites of the segment [1, £]. 

We now apply O a ^ to the state Eq. ([32]) . with a < 
X- This state can be expanded , |a7')i£A 7 |7 , /3)£-|-i j 7v- 
Using the approximate orthornormality of the states on 
the segment from 1 to £, we find that the operator in fact 
transforms a\ into a-i- Intuitively, the sum over 7, the 
state of the internal end of the £— site segment ensures 
that this operator keeps the right type of entanglement 
between the left and right side of the "cut" at £. The 
error of (|A4p scales as F(\)e~ ? , where F (\) is a function 
of x- (F(x) grows with Xi hence, to deal with a larger 
range of "energies", a larger value of £ must be used. 
This is because, as the "energy" of a state increases, it 
penetrates further into the bulk.) 

Note that A f^ does not transform high- 
entanglement energy states correctly, that is states 
such as |a/3)uv with a > x- This is because |«i7) and 
1 07') arc orthonormal only if they are both low energy 
states. This limitation is unavoidable: it is not possible 
to find a perfect representation for an operator, such as 
Q, in terms of just the £ sites near each end. One can 
add an extra particle somewhere outside of the reach 
of these sites, changing the value of Q but not of an 
observable on the ends of the chain. The physical energy 
of this state may not be much greater than the gap. 
However, not being able to describe states like this is 
not a problem when one is studying the ground state 
of the system: its entanglement energy is large, which 
means that it contributes negligibly to the value of any 
observable in the ground state. 



Appendix B: Factoring antiunitary operators 

In the analysis of time-reversal symmetry, we defined 
a parameter n by factoring T into two operators U A and 
U B , acting near the two opposite edges of the segment. 
To determine how chains add to one another, it is nec- 
essary to know the commutation and anticommutation 
properties of these operators. We ignored a small detail, 
however: since T is antiunitary, it cannot be factored 
cither as the product of two unitary or two antiunitary 
operators. One solution is just to explicitly write how T 
transforms the basis states as we did in Eqs. (fTTj) and 
(fT2"| . This becomes cumbersome after a while however, 
and later in the text we have treated U A and U B as uni- 
tary operators in Hilbcrt space, without keeping explicit 
track of their indices. Here we will explain the meaning 
of this. 



We will first discuss the bosonic case. Eq. (|11[) gives 
the action of T only on basis states. Taking a superpo- 
sition gives a factorization of T that is correct for any 
state: 

T = U A U B K, (Bl) 

where XJ a and Ub are unitary operators at the two ends 
and K is defined by 

Kj2ac<e\aP)=J2a* a0 \aP), (B2) 

a/3 a/3 

where a a p are arbitrary coefficients. Wc can now de- 
fine k by (U A K) 2 = e iK t. Thus, U A K is an antiunitary 
symmetry squaring to —1 in the nontrivial phase as in 
Kramers' theorem. Note that this equation is equivalent 
to the definition given above, Eq. (|13p . This is because 
KU A K~ 1 = U A * when the matrices are represented in 
the basis \aj3). (Note that the complex conjugate of a 
matrix depends on the basis being used, unlike the ad- 
joint.) 

We can argue physically that operators U A and U B 
satisfying Eq. (|B1| can always be found. Consider the 
ratio TK~ X between T, which is represented by complex 
conjugation in terms of the microscopic degrees of free- 
dom, and K which describes complex conjugation in the 
entanglement eigenstate basis. This operator is unitary. 
Furthermore it acts independently on the two ends: one 
may check that KO A ' B K is an operator acting on end A 
or B respectively, by expressing it in the basis of entan- 
glement cigenstates. 

The operator U A K used to define k is non-local. It 
does not commute with operators at end B, because it 
takes complex conjugates of them. However, we can 
still argue that k is well-defined: Square Eq. (|B1|) : 1 = 
U A U B KU A U B K. Since K 2 = 1, we can write this also 
as 1 = U A U B {KU A K){KU B K). KU A K is an oper- 
ator which acts on end A, therefore it commutes with 
U B . Hence 1 = [U A (KU A K)][U B (KU B K)}. Since the 
two factors are local, each must be a pure phase, hence 
(U A K) 2 = e iK l. 

Now the operators U A and U B are not uniquely de- 
fined because complex conjugation, K, is basis depen- 
dent. Changing the basis of cigenstates in which Eq. (|B2[) 
is imposed (e.g., multiplying the entanglement states by 
phase factors) changes K. This does not change topolog- 
ical properties like the value of k however: the unitary 
transformation that changes the basis can be carried out 
continuously, starting from the identity. In this process, 
k cannot change because it can only be or tt. 

For fermionic systems with /i — 0, one can decompose 
T using Eq. (|B1[) . When \x = tt, the situation is more 
complicated because the parity eigenvalue q cannot be 
associated with either one of the edges. We have to make 
sure that K still maps operators at each end of the system 
to other operators at that end. In particular, KQ A K 
must be a local operator at end A. 

This condition is satisfied if K is defined to be com- 
plex conjugation in the basis \ct(3q) provided that Q A 
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and Q are represented by either purely real or purely 
imaginary matrices in that basis. One way to satisfy 
this requirement is to first choose a basis for q = +1 
and then to construct the states in the sector from them, 
\a(3, q = — 1) = Q A \af3, q = +1). Then Q A is represented 
by u x , acting in the q = ±1 basis. In this basis, each state 
is an eigenstate of Q = <j z . Last, Q B = iQQ A = a v . 
Since Q A is real and Q B is imaginary, the two ends are 
not mixed by applying K. (If the relative phases of the 
basis states are changed, then simple complex conjuga- 
tion would mix Q A and Q B into one another.) We have 
taken the convention that Q A is real and Q B is imaginary 
above. 

Now let us show how to calculate k when two 
(/!, (f>, k) = (-7T, 0, 0) chains are combined. The factor- 
ization T = U A U B K must be carried out in a basis of 
states of the form |a/3), according to our conventions. 
One basis for the states on the two chains together is 
given by {|±)i| ± ) 2 } (where the sign represents the val- 
ues of qi, <72- (We do not explicitly write the bosonic in- 
dices a, /3.). These states map to themselves under time 
reversal. However, U A and U B cannot be the identity be- 
cause we know they must be fermionic; this is the wrong 
basis for defining K by simple complex conjugation. 

Let us transform the states to a basis in which there 
is no entanglement between the ends; we therefore use 
states that are eigenvectors of the local operators Q A = 
-1Q1Q2 and Q B = iQ B Q B , namely |?a(?b>- (The rel- 
ative minus sign between Q A and Q B ensures that the 
total parity is qAlB = 9i92-) 

To construct the basis, first find an eigenfunction of 
Q A and Q B with eigenvalues +1, (choose the phase ar- 
bitrarily): 

|+A+ B > = ^(|+)i|+>2-i|-)i|-> 2 ) (B3) 

Now generate the other basis states from this by applying 
Q A and Qf : 

\-A+B) = Qf|+A+B) = ^(|-)l|+) 2 -i|+)l|-> a ) 
\+A-B) = -igf|+A+fl> = ^(|->l|+>2+<|+>l|-> a ) 
\-A-B) = Qf\+A-B) = ^=(\+)l\ + )2+ l \-)l\-)2) 



(B4) 

The phases are just conventions in the first two defini- 
tions, and the phase in the third equation follows from 
the independence of the ends: Q A has to act on qA the 
same way no matter what the value of qs is. Now it 
is clear that T switches the fermion parity of each end 
in this basis, since changing the sign of i exchanges the 
states \qAqB) and | - qA, -gs)- 

Now we can define K to map each of these basis states 
to itself. Clearly, T = Q2 because Q2 also exchanges the 
same pairs of wave functions, or more precisely T = Q2K . 
Hence U A = -iQ A and U B = Q B . One can check that 
{U A K f = -1, so k = tt. 

Note that, in spite of all this trouble, the value of k in 
a phase with fi = 0, <f> — n does not have any physical 
significance-the four-fold degeneracy of the spectrum is 
already explained by the fact that 4> — tt. The reason 
T changes fermion parity at the ends is that the two 
ends can only be disentangled by a change of basis that 
includes complex phases. 



Appendix C: Parity of the ground state on a 
periodic chain 

The parity of the ground state on a periodic chain is 
given by e l ( eL +^) ^ This follows from a fact in Sec. [V] 
Consider two subsegments of the chain, one ending at X 
and the other starting at X + 1. The ground state wave 
function is an eigenfunction of the following: 

e~^Q B (X)Q A (X + lM = \^). (CI) 

When fi — 7r, this relation describes the correlations be- 
tween the Majorana degrees of freedom in adjacent seg- 
ments of the chain. 

We now suppose the periodic chain has length L and 
break it at two places, between V and L' + 1 and 
between L and 1. The total fermion parity of the 
ground state is the product of the parity on the two seg- 
ments, (e l ^ L '-^Q A (l)Q B (L'))(e^ L - L '^-^Q A (L' + 
1)Q B {L)). Rearranging and using Eq. (|C1J) . the ground 
state parity comes out as e l ( eL +^) . The extra minus sign 
when fi = 7r comes from anticommuting the Q operators. 
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